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The laser induced fluorescence (LIF) spectra yFE + ~ & 3 II(iJ J ) — > X 1 E + of KCs dimer were 
recorded in near infrared region by Fourier Transform Spectrometer with a resolution of 0.03 cm -1 . 
Overall more than 200 collisionally enhanced LIF spectra were rotationally assigned to 39 K 133 Cs and 
41 K 133 Cs isotopomers yielding with the uncertainty of 0.003-0.01 cm -1 more than 3400 rovibronic 
term values of the strongly mixed singlet A 1 E + and triplet & 3 II states. Experimental data massive 
starts from the lowest vibrational level va = of the singlet and nonuniformly cover the energy 
range E J £ [10040,13250] cm" 1 with rotational quantum numbers J £ [7,225]. Besides of the 
dominating regular A 1 E + ~ t> 3 IIn=o interactions the weak and local heterogenous J 4 1 E + ~ & 3 IIn = i 
perturbations have been discovered and analyzed. Coupled-channel deperturbation analysis of the 
experimental 39 K 133 Cs e-parity termvalues of the A 1 E + ~ & 3 IIn = o,i,2 complex was accomplished 
in the framework of the phenomenological 4x4 Hamiltonian accounting implicitly for regular in- 
teractions with the remote 1 I1 and 3 E + states manifold. The resulting diabatic potential energy 
curves of the interacting states and relevant spin-orbit coupling matrix elements defined analytically 
by Expanded Morse Oscillators model reproduce 95% of experimental data field of the 39 K 133 Cs 
isotopomer with a standard deviation of 0.004 cm" 1 which is consistent with the uncertainty of the 
experiment. Reliability of the derived parameters was additionally confirmed by a good agreement 
between the predicted and experimental termvalues of 41 K 133 Cs isotopomer. Calculated relative in- 
tensity distributions in the A ~ b — > X LIF progressions are also consistent with their experimental 
counterparts. Finally, the deperturbation model was applied for a simulation of pump-dump optical 
cycle a 3 E + — > J 4 1 E + ~ fe 3 II — > X 1 E + proposed for transformation of ultracold colliding K+Cs pairs 
to their ground molecular state vx = 0; Jx — 0. 

PACS numbers: 



I. INTRODUCTION 

In ultracold (below 1 mK) gases of polar alkali 
diatomic molecules the electric dipole-dipole inter- 
action possesses long-range and anisotropic charac- 
ter and cause particular attraction since they are 
promising for a variety of novel applications 
Among them there are proposals of quantum phase 
transitions 0, quantum information devices [||, 
coherent control of chemical reactions [|[ ■ Further 
success in producing ensembles of ultracold quan- 
tum gases of diatomic molecules relies on the accu- 
rate knowledge of the molecular rovibronic struc- 
ture and transition probabilities gained from the 
high resolution spectroscopy and state of art ab 
initio calculations. 

As far as polar alkali diatomic molecules arc con- 
sidered, experiment based information on the po- 
tential energy curves of the lowest excited A 1 E + 
and 6 3 n states mixed by singlet-triplet spin-orbit 
(SO) interaction is of great importance due to the 



possibility of their usage as intermediate states for 
transferring the vibrationally excited molecules ob- 
tained from cold colliding atoms via photoassocia- 
tion or Feshbach resonances into the absolute rovi- 
bronic ground states X 1 E + with vx = 0, Jx = 
0, Q . In the sequence of heteronuclear alkali di- 
atomics KCs is among promising species for pro- 
ducing ultracold polar quantum gases. Due to 
proximity of K(4p) and Rb(5p) energies, the en- 
ergy of low-lying electronic states of KCs is sim- 
ilar to the ones in the RbCs molecule, however, 
KCs possess larger permanent electric dipole mo- 
ment. The RbCs molecule has been already stud- 
ied in ultracold conditions 0, H H- In particu- 
lar, the recent quantum dynamics simulation of 
the mixed A l Y> + and 6 3 IIq states of RbCs in a 
time-dependent wave-packet approach Q demon- 
strated that the pump - dump picosecond pulsed 
laser scheme could be efficient to form ultracold 
heteronuclear diatomic molecules in highly bonded 
ground X 1 Y> + state. The highly accurate em- 



pirical potentials for the ground singlet A 1 Yj + 
and triplet a 3 £ + states converging to the lowest 
K(4S)+Cs(6S) asymptote (see Fig. QJ as well as 
predicted scattering lengths and Fcshbach reso- 
nances which are required to simulate cold collision 
processes are presented in Ref's |Tl|] . However, 
till now there is no empirical information about 
the excited states of KCs while ab initio potential 
energy curves (PECs) have been calculated for a 
wide range of internuclear distance in the frame- 
work of pure Hund's (a) QJ, [H, 03 and (c) [H 
coupling cases. The relevant permanent and tran- 
sition dipole moment functions are also available 
from Ref's [HQ. 

Particular complicated objects for the theoreti- 
cal interpretation are alkali diatomics containing a 
heavy Rb or Cs atoms since the relevant singlet- 
triplet SO coupling matrix element £jjg> between 
A 1 Yi + and b 3 Hn=o states is comparable with vibra- 
tional spacing of the interacting states. Therefore 
y4. 1 S + and 6 3 rin=o states are fully mixed (later de- 
noted as the A ~ b complex) in both adiabatic 
and diabatic basis set representation [l6[. First 
studies of such kind of system have been reported 
for the Rb-containing alkali diatomics, in which 
apparent disorder in vibrational spacing was ob- 
served for Rb-? fl7l . see [lU for more extensive 
studies, NaRb [3 and RbCs [H. The crucial is- 
sue however remained unsolved, namely, while the 
experimental term values were obtained with high 
precision of about 0.003-0.01 cm -1 , the resulting 
depcrturbed parameters reproduced experimental 
data with much poorer accuracy of 0.05-0.25 cm" 1 
[3 O Furthermore, even the vibrational 

numbering of the "dark" triplet b 3 H state was still 
remained questionable. 

In recent studies on the A ~ b complex of NaRb 
[2~i| and NaCs [22| more comprehensive depertur- 
bation method has been developed in the frame- 
work of the inverted channel-coupling approach by 
means of the 4x4 Hamiltonian for the A 1 Yi + and 
6 3 rio=o,i,2 substates constructed on Hund's cou- 
pling case (a) basis functions. Phenomenological 
inclusion of the splitting of the 6 3 nn=o,i.2 sub- 
states and of indirect A 1 Y> + and 6 3 nn=i states 
coupling allowed accounting for the regular per- 
turbation. The elaborated model reproduced the 
experimental term values obtained by the high res- 
olution Fourier transform spectroscopy (FTS) of 
the laser induced fluorescence (LIF) from/to the 
A ~ b complex with a standard deviation of 0.006- 
0.012 cm" 1 , which was consistent with the uncer- 
tainty of the FTS experiment. The paper [22] 
presents the most challenging though attractive, 



both for model's testing and applications, case of 
A ~ b complex study of a Cs-atom containing al- 
kali diatomics where SO interaction is of largest 
value mainly determined by £ so ~ 185 cm" 1 for 
6 2 P(Cs). This successful approach opens the path 
to study of the A ~ b complex of more heavier 
Cs-containing heteronuclear alkali diatomic such 
as RbCs and KCs. While the RbCs molecule is 



currently under investigation 23j the present pa- 
per deals with the KCs molecule. 

It should be noted that the present experimental 
and deperturbation studies of the KCs A ~ b com- 
plex are greatly facilitated by the couple of facts. 
First, the KCs molecule was recently involved in 
detailed FTS LIF investigation [l(| which resulted 
in high accuracy empirical point-wise PEC for the 
ground X 1 T, + state covering a range vx — to 
97 and Jx = 12 to 209. Moreover, a majority of 
the data for PEC construction in [10[ was obtained 
directly from the (A ~ b) — > X LIF spectra con- 
taining the information on rovibronic energies of 
the A ~ b complex. Second, the accurate ab initio 
spin-orbit coupling matrix elements between the 
lowest electronic states of KCs have been obtained 
recently in Ref. [lij] as explicit functions of inter- 
nuclear distance. 

The goal of the present work is to perform high- 
resolution spectroscopy studies of the A ~ b com- 
plex of the KCs molecule and to apply the ad- 
equate data processing in order to get empirical 
molecular structure parameters, namely the deper- 
turbed A 1 S+ and & 3 n diabatic potential energy 
curves along with the related SO coupling matrix 
elements functions on internuclear distance R. 

In summary, Sections II and III present the ex- 
perimental techniques and spectra analysis, Sec- 
tion IV presents the deperturbation model and 
computational details, Section V includes obtained 
fitted parameters, as well as interpretation of ex- 
perimental data. The paper ends with short con- 
clusion (Section VI). 



II. EXPERIMENT 

In the experiment the rovibronic levels of the 
KCs A ~ b complex were directly exited by diode 
lasers and subsequent (A ~ b) — > X LIF spectra 
were recorded by Fourier transform spectrometer. 
The experimental setup was the same as used for 
KCs ground state studies [Tcj and, hence, we will 
mention here only some key points of the exper- 
iment. KCs molecules were produced in a linear 
stainless steel heat pipe. The heat pipe was filled 
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FIG. 1: (Color online) Scheme of electronic terms 
based on ab initio results of Ref. [l^ ] 

with 10 g of K (natural isotope mixture) and 7 
g of Cs (both metals in ampules from Alfa Ae- 
sar). The typical operating pressure of Ar buffer 
gas in the heat pipe was 3-5 mbar. On the one 
hand Ar gas prevented condensation of the metal 
vapor on the heat pipe windows, and, on the other 
hand, collisions of optically excited KCs molecules 
in particular rovibronic levels v' , J' with Ar atoms 
ensured population of great number of neighboring 
rotational levels and observation of the respective 
transitions to the ground state. During the ex- 
periments the heat pipe was kept at 270 — 290°C 
temperature by a Carbolitc furnace. 

In the experiment the laser beam was sent into 
the heat pipe through a pierced mirror. Backward 
LIF was collected by the same mirror and focused, 
by two lenses, on the input aperture of the spec- 
trometer (Bruker IFS 125HR). For LIF detection 
we used an InGaAs diode operated at room tem- 
perature. The spectral sensitivity of the detec- 



tor gradually diminishes from about 6500 cm" 1 
toward higher frequencies reaching 60% of maxi- 
mal sensitivity at ~ 10000 cm" 1 . Then sensitiv- 
ity drops much faster, reaching 10-15% of maximal 
value at 11000 cm" 1 . The resolution of the spec- 
trometer was typically set to 0.03 cm . 

The (A ~ b) <— X transitions were excited by 
tunable diode lasers. Four single mode laser diodes 
with central wavclcnghts 850 nm (LD850/100 from 
Toptica Photonics), 980 nm (L980P200I from 
Thorlabs), 1020 nm (LD-1020-0400 from Toptica 
Photonics) and 1060 nm (L1060P100J from Thor- 
labs) were mounted in homemade external cav- 
ity resonators (Littrow configuration) with a grat- 
ing serving as a feedback source. The respec- 
tive frequency tuning ranges were 11560 - 11930 
cm" 1 , 10209 - 10515 cm" 1 , 9700-9860 cm" 1 , and 
9360-9510 cm" 1 . Fine tuning of the grating was 
achieved by a piezoelectric actuator. Temperature 
and current stabilization was ensured by Thorlabs 
controllers. The power of the lasers at the entrance 
of the heat pipe varied from 15 to about 50 mW 
depending on particular laser diode and laser cur- 
rent. The general approach to measurements was 
the following. The laser frequency was tuned until 
the LIF signal monitored at reduced resolution in 
the "Preview Mode" of the spectrometer gets its 
maximal value and then fixed during recording of 
the spectrum. In order to ensure sufficient signal- 
to-noise ratio for the lines of medium strength, the 
number of scans for each recorded spectrum varied 
from 20 to 40; averaging over a number of repeated 
measurements was applied in some cases. In some 
cases the strong scattered laser light from the in- 
put window of the heat pipe was eliminated by 
long-pass edge filters. 

III. ANALYSIS OF THE SPECTRA 

Recorded KCs LIF spectra, see Figs. 2-5, typi- 
cally consisted of one or two strong doublet pro- 
gressions and several weaker ones. 

The assignment of the LIF progressions in the 
(A ~ b) — ► X spectra was straightforward 
from measured vibrational and rotational spacings 
thanks to accurate ground state PEC [lj|. The 
term values of the upper rovibronic levels of the 
A ~ b complex giving rise to (A ~ b) — > X LIF 
were obtained by adding the corresponding ground 
state level energy to a transition wave number. 
The uncertainty of the line positions is estimated 
to be 0.1 of the resolution, or 0.003 cm" 1 . For 
lines with signal-to-noise ratio (SNR) less than 3 
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FIG. 2: 39 K 133 Cs (A ~ b) -> X LIF spectrum 
recorded by InGaAs detector at excitation with 850nm 
laser diode (laser frequency 11608.623 cm -1 ). The 
strongest doublet P63 , Rai progression 5 < vx < 78 
originates from the level J' = 62 with energy E — 
12164.984 cm -1 excited by the transition (i>*,62) *— 
(6,63). Another weaker KCs LIF progression excited 
in («*,133) <— (6,134) transition is not marked. A 
mesh of weaker lines around 9500 cm -1 belongs to K2 
LIF. 



the uncertainty increases and was assumed to be 
about 0.01 cm -1 for lines with SNR of about 2. 
The term values of the upper state levels were de- 
termined as average value from at least four transi- 
tions. Hence, the main uncertainty of the term val- 
ues could be caused by the shift (if any) of the laser 
frequency from the center of the Doppler broad- 
ened (A ~ b) <— X absorption line. When LIF is 
registered, as in our case, along the direction of the 
laser beam, this shift is transferred to the respec- 
tive shift of the LIF line positions. The Doppler 
broadening for KCs at the exploited spectral region 
is about 0.012 cm -1 and we estimate maximal ex- 
perimental uncertainty of the A ~ b complex term 
values as 0.01 cm" 1 . 

Typically, around the strong LIF lines the satel- 
lite lines could be recorded due to collision - in- 
duced distribution of the population of the directly 
excited rovibronic level over neighboring rotational 
levels and even vibrational levels. The assignment 
of these lines was based on the high accuracy of the 
ground state PEC and it was always checked that 
the spacing of satellite doublet lines (vr — Vp) co- 
incides with calculated energy differences between 



the respective ground state rotational levels within 
accuracy of some miliwavenumbers. 

Figure 2 represents an example of the recorded 
LIF spectrum when the (A ~ b) <— X transition 
is excited with the 850 nm laser diode tuned to 
11608.623 cm -1 . The recorded spectral line inten- 
sities for frequencies above 10000 cm -1 are sub- 
stentially diminished because of dropping of the 
spectral sensitivity of the InGaAs detector. The 
LIF intensity distribution in such a long progres- 
sion as depicted in Fig. 2, with strong last max- 
imum at high vx, is characteristic for highly ex- 
cited vibrational levels va of the A ~ b complex. 
Note that, in general, the LIF spectra recorded 
with 850 nm diode excitation usually contained a 
very strong A 1 T, U — ► AT 1 S ff band of K 2 , and only 
by careful selection of the laser frequencies it was 
possible to eliminate the K2 fluorescence, at the 
same time having strong enough KCs (A ~ b) — > X 
LIF progressions. 

The usage of 1020 nm and 1050 nm laser diodes 
allowed us to reach low lying A - state rovibronic 
levels of KCs. Rather strong LIF signal made it 
possible in a number of cases to detect spectra 
without filtering of scattered laser light. Thus, full 
LIF progressions could be recorded, which was im- 
portant for establishing the relative intensity dis- 
tribution in the progressions, which is helpful for 
vibrational assignment of the upper state. The ex- 
ample of such a spectrum is given in Fig. 3. This 
spectrum contains five assigned KCs (A ~ b) — > X 
progressions, two of them with J' = 50 and 109 
are originating from the lowest va level with pre- 
dominant singlet character. These two progres- 
sions show unusual intensity distribution. Along 
with one strong bell - shape maximum in inten- 
sity distribution around 9500 cm -1 at low vx, as 
is typical for progressions from upper state v' = 
vibrational level, we have found for these progres- 
sions transitions to higher vx resulting in addi- 
tional intensity maximum, however substantially 
weaker than the main one. The small fragment of 
the spectrum with these additional transitions is 
given in Fig. 3b. Besides, the spectrum in Fig. 3 
contains a progression originating from the level 
with energy E' = 10227.038 cm -1 and J' = 44 
which is of predominantly triplet character as will 
be shown below, see Section V. A short fragment 
of this progression can be seen in Fig. 3b where 
the respective P, i?-doublets are marked above the 
spectrum. In general, the laser most often exited 
the levels with a dominant singlet character, and 
only in some cases a direct excitation of rovibronic 
(A ~ &)-state levels with a dominant triplet char- 
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FIG. 3: (a) The 39 K 133 Cs (A ~ b) -» X LIF spec- 
trum recorded at laser frequency 9800.135 cm" . The 
following (j4 ~ 6) <— X; (v*,J r ) <— (vx,Jx) excita- 
tion transitions have been assigned: («*, 50) <— (4, 51); 
(«*,109) <- (3,110); 0*,57) <- (7,56); (u*,123) <- 
(16,122); (i>*,44) <- (5,43). Weak Cs 2 progressions 
have been assigned below 8500 cm" 1 . An indexes for 
Pj x , Rj x lines denote the respective ground state ro- 
tational levels, (b) The zoomed fragment of the spec- 
trum; doublets originating from the lowest excited vi- 
brational level posessing dominant singlet character 
with J' = 50 and J' = 109 are marked with black 
and red vertical bars below the spectrum. The bars 
above the spectrum mark transitions from the level 
(u* , J' = 44) with dominant triplet character. Num- 
bers at bars denote vx- 
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FIG. 4: (a) The 39 K 133 Cs {A ~ b) -> X LIF 
spectrum recorded at excitation frequency 9743.083 
cm— 1. The scattered laser light and LIF at fre- 
quencies higher than about 9500 cm— 1 are cut 
by a long-pass edge filter. The following excita- 
tion transitions for 39 K 133 Cs have been assigned: 
(u*,49) <- (6, 50); (v* , 81) <- (17,82); («*, 25) <- 
(13,26); (v*, 135) <- (12, 134); («*, 150) <- (11,149). 
(b) Zoomed part of the spectrum around vx = 20 for 
the strongest progression with J' = 49. Three groups 
(II-IV) of satellite P, R branches shifted from the main 
rotational relaxation pattern (group I) are clearly seen. 
The inset presents term values of the respective rovi- 
bronic levels as dependent of J(J+1). Directly excited 
level J' = 49 ir marked. 



acter took place. 

Fig. 4 and Fig. 5 demonstrate how fruitful is 
the collision induced rotational relaxation for sys- 
tematic study of the A ~ b complex level struc- 
ture. The spectrum in Fig. 4a contains one very 
strong progression originating from the level E' = 
10255.413 cm" 1 , J' = 49 as well as five weaker 
ones. A large number of satellite lines was as- 



signed around the strong P^,R^ doublet lines, 
especially at transitions to vx = H- The most in- 
teresting feature consists in appearance of the ad- 
ditional three wide groups of P, R satellite lines, as 
seen in Fig. 4b. These groups are shifted from the 
main P50 , i?48 doublet and represent a situation 
when an initial population of the optically excited 
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FIG. 5: (Color online) (a) The LIF spectrum recorded 
at excitation frequency 9826.348 cm - . The strongest 
progression originates from the upper level excited in 
transition (v* , 158) <— (4, 157). (b) zoomed part of the 
spectrum around vx = 12; Red and blue bars denote 
assigned P and R satellite lines respectively. Long bars 
mark transitions from upper state levels with dominant 
triplet b 3 IIi character. Numbers in the spectrum de- 
note corresponding ground state rotational quantum 
numbers for observed transitions. 



level with J' = 49 is collisionally distributed within 
two state ( J 4 1 E + and fr 3 IIo) rovibronic level mani- 
folds of the A ~ b complex. As a result, more than 
130 term values were obtained from this spectrum, 
see the inset in Fig. 4b. Note that these satellite 
groups are the strongest in the range 8700 cm" 1 
- 9000 cm -1 and substantially weaker at higher 
frequency (lower vx) range. 

In several spectra clear evidences of local per- 
turbations were observed as illustrated in Fig. 5. 
The spectrum given in Fig. 5a is recorded at ex- 
citation frequency 9826.348 cm -1 ; it consists of 



FIG. 6: (Color online) Experimental term values of the 
(A X E+ ~ 6 3 n) complex plotted versus J'(J' + 1). 



one strong KCs LIF progression from the level 
(E' = 10861.046 cm" 1 , J' = 158) and several weak 
KCs and K 2 LIF progressions, hardly seen in this 
scale. Zoomed part of the strong progression in 
the range of vx = 12 is given in Fig. 5b. Along 
with an ordinary satellite line pattern around the 
central P, R lines with gradually changing spac- 
ing between the lines, the local perturbation can 
be clearly recognized from the two gaps ("win- 
dows") in the both P,R branches, see P159/P158 
and i?i57/i?i56- Moreover, we have assigned "ex- 
tra" lines with the same rotational quantum num- 
bers, namely P 158 , P 159 , P' 156 and R' 157 , see Fig. 5b. 
As it will be shown later these additional lines take 
their origin from the levels with dominant 6 3 IIi 
state character and they can be seen here due to 
effective collisional transfer of population from the 
directly excited predominantly singlet A state level 
with J' = 158, since it is very close to the pertur- 
bation center. 

Overall such local perturbations accompanied 
with extra lines were observed in three spectra. In 
several cases the local perturbations were observed 
as irregularities in line sequences without "extra" 
lines. 

Summarizing, during this study more than 3400 
term values of the A ~ b complex rovibronic lev- 
els were obtained. These term values are given in 
Fig. [5] as dependent on J' (J' + 1). Empty circles 
seen in Fig. [5] show distribution of 31 term value 
determined for the 41 K 133 Cs levels. The data for 
41 K 133 Cs are given also in Table HI 



6 



TABLE I: The experimental rovibronic termvalues 
E expt (in cm -1 ) of the of 41 K 133 Cs isotopomer assigned 
to the A 1 E+ ~ 6 3 n complex. A = E expt - E cc is 
the difference between experimental and predicted en- 
ergies. Pi (in %) are fractional partition of the levels, 
see Section V. 



J' E expt A P A P b0 P bl 

59 10205.767 0.003 44.9 55.1 0.0 

75 10250.792 -0.001 85.2 12.7 2.1 

58 10424.388 0.001 65.8 34.2 0.0 

52 10547.339 -0.008 66.3 33.5 0.2 

53 10549.850 -0.008 66.8 32.9 0.3 

54 10552.405 -0.007 67.2 32.3 0.5 

55 10555.010 0.000 67.5 31.6 0.9 

56 10557.639 -0.014 67.5 30.7 1.8 

62 10555.967 -0.003 51.1 48.8 0.1 

63 10559.120 0.001 50.4 49.5 0.1 

64 10562.317 -0.007 49.7 50.2 0.1 

65 10565.579 -0.007 49.0 50.9 0.1 

66 10568.900 -0.006 48.3 51.6 0.1 

67 10572.278 -0.005 47.6 52.4 0.1 

68 10575.711 -0.006 46.8 53.1 0.1 

69 10579.199 -0.010 46.1 53.8 0.1 

70 10582.754 -0.007 45.3 54.6 0.1 

71 10586.355 -0.014 44.6 55.3 0.1 

72 10590.031 -0.006 43.9 56.1 0.1 
69 10595.554 -0.006 75.6 24.3 0.1 
62 10713.205 -0.006 61.5 38.4 0.1 
25 10496.360 -0.011 55.7 44.3 0.0 

73 10608.598 0.005 77.4 22.6 0.1 
34 10575.252 0.002 83.9 16.0 0.0 
87 10294.767 0.004 86.3 13.6 0.1 
20 10352.114 0.004 61.3 38.6 0.0 
119 10662.753 -0.001 77.3 22.5 0.1 
143 10801.142 -0.001 70.9 27.4 1.8 
134 10968.075 0.005 80.2 19.5 0.3 
49 12053.926 -0.001 82.7 17.3 0.0 
71 12511.108 0.009 81.7 18.3 0.0 



IV. DEPERTURBATION ANALYSIS 
A. Modeling Hamiltonian 

The non-adiabatic rovibronic wave-function 

*j = Ei<t>m (* e [A^, & 3 n fi=0 ,i, 2 ]) cor- 

responding to the j-th rovibronic level of the 
A 1 E + ~ & 3 I1 complex of KCs molecule with 
the fixed rotational quantum number J and e- 
symmctry was approximated by the linear com- 
bination of the symmetrized electronic-rotational 
wavefunctions ifi belonging to a pure Hund's cou- 
pling case (a) The expansion coefficients 
<)>i (r) are the r-dependent fractional components of 



the non-adiabatic vibrational wavefunction <&j = 
\\<t>Ai 4>bo, 061 j 0&2 1 1 which are determined by the 
bound solution of the close-coupled (CC) radial 
equations [M HI HH : 

{- 1 ^ +Y{r ^' J) - 1E ? C )^ {r)=0 (1) 

with the conventional boundary </>i(0) = 4>i(oo) = 
and normalization Pj = 1 conditions, where 
Pi = (^i\<f>i) is the fractional partition of the j- 
th level. Here I is the identity matrix, E^ c is the 
total non-adiabatic energy of rovibronic level of the 
complex while V is the symmetric 4x4 matrix of 
potential energy consists of the diagonal 

Vi s+ = U A + B[X + 2] (2) 

^n = U b0 + B[X + 2] 

V» Ul = U bl + B[X + 2] 

Vsn 2 = U b2 + B[X-2] 

and non- vanishing off-diagonal 

V^ + -m a = -Vk%o (3) 

V3 Uo _ mi = -BV2X 

ra ni -»n a = -B^2(X-2) 

^s+- 3 ni = —BCAbW2X 

matrix elements explicitly depending on reduced 
mass (i and rotational quantum number J as a 
parameter: 

Hereafter all electronic parameters of the model 
are assumed to be mass-invariant. 

Here Ua (r), Ubd(r) are the diabatic PECs of 
the singlet and triplet 6 3 nn = o ! i.2 sub-states 

while £%° bo (r) is the relevant spin-orbit coupling 
matrix element. Besides of strong off-diagonal ho- 
mogenous A 1 Y, + ~ 6 3 nji=o interactions the de- 
perturbation model also explicitly consider het- 
erogenous spin-rotational interaction between all 
H = 0,1,2 components of the triplet 6 3 n state 
[l6| . Moreover, the non-equidistant SO splitting 
functions ^41° ^ between the 6 3 IIn = o.i.2 sub- 
states 

U bl = U bQ + Al°; U b2 = U bl + A s ° (4) 

and the indirect coupling A 1 H + ~ b Hn=i param- 
eter t^Abi were introduced in the Hamiltonian ([2]); 
([5]) in order to account implicitly for regular per- 
turbations caused by the remote 1 n and 3 S + states 
manifold (first of all the _B 1 n and c 3 S + states con- 
verging to the second dissociation limit, see Fig.[T|). 



7 



B. Computational details/Fitting procedure 

The empirical PECs of the interacting A 1 E + and 
6 3 n states Ua{t), Ubo(r) as well as both diago- 
nal A±(r) and the off-diagonal £,%° b0 (r) spin-orbit 
functions were represented analytically by the Ex- 
panded Morse Oscillator (EMO) function [H[ 




(5) 



ref 



rP + r p ref 



converging to the appropriate atomic limit. In 
particular, for all SO functions T*° = - S*° 
where fg> s = [E 6 2p 3/2 - E 6 2 Pl/2 }/3 is the spin- 
orbit constant of the Cs atom in the 6 2 P state 
[13, HI- For diabatic PEC U A of the singlet state 
Tf = T dls - while = T dls - &° - D b e ° 
for the Ubo sub-state of the triplet. Here Tdis — 
+ E 6 2 P — E 6 2 S is the energy of center of gravity 
of the Cs(6 2 P) doublet (without hfs splitting) [27| 
referred to the minimum of the ground A-state. 
The Df = 4069.3 cm -1 value is taken from Ref. 

The initial EMO parameters (£> e , r e and a,i) of 
the relevant PECs were estimated by averaging of 
the independent ab initio results borrowed from 
Ref 's [H, 03 while the required SO functions 
were constructed using the recent quasi-relativistic 
estimates of Ref. Off-diagonal SO function 

could be empirically determined only in the close 
vicinity of the crossing point of the diabatic PECs 
of the interacting A 1 Yj + and 6 3 I±o states. There- 
fore, the required SO-EMO functions were con- 
strained by fixing of the respective a,(i > 0) pa- 
rameters based on the relevant ab initio points [14| , 
so only Df, r s ° and af° SO-EMO parameters were 
rest to be variable in the fit. 

The refined parameters of the PECs and SO 
functions combined with the r-independent £ A °bi = 
const parameter were determined iteratively dur- 
ing the weighted nonlinear least-squared fitting 
(NLSF) procedure procedure 



xt xpt 



N, 

E 

0=1 



Wj (Et 



cc 



riexpt 
3 



N, 



expt 



(6) 



where Wj = 1 /<7 2 is the weight of each level, E^ xp 
its experimental termvalue and o~^ xpt = 0.003—0.01 
cm -1 its uncertainty. Here, N expt is the number 
of experimental termvalues involved while M p is 



the total number of adjusted parameters of the 
model. Only experimental rovibronic termvalues 
gexpt o j mos t abundant 39 K 133 Cs isotopomer 
were included in the fitting procedure. The min- 
imum of the functional (J5]) was searched by the 
modified Levenberg-Marquardt algorithm [3l| re- 
alized by MINPACK software [H. 

The analytical mapping procedure based on re- 
placement of the conventional radial coordinate (in 
A) by the reduced radial variable y{r\ f = 5.2, /3 = 
5) = [1 + (f/r) 13 }" 1 was used to transform the 
initial CC equations |T|) into completely equiva- 
lent form. Then, the modified CC equations given 
explicitly in Ref. 13 3| were solved on the inter- 
val r e [2.5, 15.5] Aby the finite-difference (FD) 
boundary value method [34j with the fixed number 
of uniformed grid points N . The central 5-points 
FD approximation (FD5) of the kinetic energy 
term was employed in Eq. |T]) . The ordinary eigen- 
value and cigenfunction problem of the resulting 
symmetric band matrix was iteratively solved by 
the implicitly restarted Lanczos method realized 
in ARPACK software in the shift-inverted spectral 
transformation mode (35| . The energy error cor- 
rection for the FD5 method recently invented in 
the explicit integral form [22| was used to extrapo- 
late to infinite number of grid points the eigenvalue 
obtained at the fixed number of integration points. 
The optimized procedure allowed us to attain the 
absolute accuracy of the calculated energies E^ c 
about of 0.001 cm" 1 using only N = 1500 - 2000 
grid points. 



V. RESULTS AND DISCUSSION 



A. Interatomic potentials and spin-orbit 
functions 



The resulting EMO parameters of the deper- 
turbed (diabatic) PECs of the and 6 3 n 

states are presented Table |H] while the respec- 
tive empirical spin-orbit functions are given in Ta- 
ble IIIII The parameters reproduce 95% of ex- 
perimental data field of the 39 K 133 Cs isotopomer 
with a standard deviation (SD) of 0.004 cm -1 in 
the framework of the phenomcnological model de- 
fined by Eq. (0) and Eq. ©. The SD value is 
well-consistent with the estimated uncertainty of 
0.003 — 0.01 cm' 1 of the present experiment. The 
EMO parameters defined by Eq. ([5]) are given in 
the EPAPS in double precision format [3g|. For 
convenience the fitted PECs and SO functions 
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TABLE II: The resulting EMO parameters of the dia- 
batic potential energy curves of the deperturbed A 1 £ + 
and 6 3 IIo states. Tdis, 2) e and T e in cm -1 , r re / and r e 
in A, m in A , p is the dimensionless. Tdis, r re f and 
p parameters were fixed in the fit. 



Parameter 




b d n 




15616.95 


p 


3 


4 


• ref 


5.0 


4.2 


B e 


5567.546 


6599.300 


r e 


4.981380 


4.179865 


a 


0.44685377 


0.56383223 


(ii 


0.01153475 


0.11731134 


11-2 


0.01224621 


0.10163399 


«:i 


0.12946458 


-0.04095288 


(1-1 


0.16407832 


-0.37280647 


tin 


0.28366171 


0.11826773 


(16 


-0.27382905 


2.22890307 


7 


-0.61363449 


0.57576966 


08 


1.05182596 


-5.39299688 


(>!> 


0.87045466 


-3.16111671 


aio 


-1.00928649 


8.86174318 


<m 


-0.72414478 


4.22763021 


(112 


0.00055024 


-7.36542364 


dl3 




-0.00010791 


ai4 




0.12543766 



arc represented in the EPAPS by the appropri- 
ate point-wise form as well. Furthermore. EPAPS 
tables contain both experimental and reproduced 
robvibronic termvalues of the complex as well as 
their residuals and fractional partitions. 

As can be seen from Fig. [7J the singlet A and 
triplet b state PECs are intersecting in the vicinity 
of the equilibrium distance r e of the A-state. As 
a result, all levels of the singlet state are strongly 
perturbed and the most significant perturbations 
take place for low vibrational levels va of the A- 
state. The diabatic (deperturbed) t^o = 18 vibra- 
tional level of the triplet state is found to be the 
first one which reach by energy the ground diabatic 
va = level of the singlet at low J'- values. How- 
ever high rotational levels of the va = basically 
interact with lower vibrational levels (vbo < 18) of 
the 6-triplet since r\ < (see Table ITV|) . 

Fig. [7] and [S] demonstrate a good agreement of 
the derived PECs and SO matrix elements with 
the ab initio quasi-relativistic results [H, [H, [3] 
both inside and outside the experimental data re- 
gion r € [3.1, 7.8] A. In particular, Table HVl clearly 
shows that the empirical molecular constants T e , 
u e and r e coincide with their ab initio counterparts 



TABLE III: The resulting EMO parameters of both 
diagonal AJ" and off-diagonal spin-orbit coupling 
functions. *& B , T) s e ° and T£° in cm -1 , f r ref and r s ° 
in A, a,i in A , T p = 1 is the dimensionless. * denotes 
the fixed parameter. The empirical r-independent cou- 
pling parameter £a&i = 0.04935. 



Parameter 


SA60 


At° A s + ° 


KCs 




184.68 


V f 

1 re j 


5.3 


5.5 


©r 


102.998 


106.419 101.539 




5.054338 


5.442732 5.576276 


a 


0.31895 


0.42157 0.37735 




0.40997 


1.02228 




0.49244 


2.23836 






-1.71726 






-9.84181 


rpSO 

* e 


81.68 


78.26 83.14 



TABLE IV: Comparison of the basic spectroscopic con- 
stants available for the KCs A 1 E + and b 3 Il states. The 
T e and co e values are given in cm -1 while r e in A. T e 
values are reffered to the minimum of the X state. The 
harmonic frequency for the present PECs U(r) is cal- 
culated as LO e = h^/k/fi [lfl, where k = d 2 U/dr 2 \ r=r(1 
is the strength constant. 



State T e r e u> e T e r e ui e T e r e uj e 

present [15] 
A L T, + 10049 4.98 49.8 
b 3 IIo 8833 4.18 68.4 8717 4.17 70.3 
& 3 IIi 8938 4.19 68.3 8856 4.16 71.6 

b 3 Yl 2 9044 4.20 68.2 8981 4.15 71.8 

[12] M [14] 

IFW 10107 4.93 51.8 9947 4.95 51.8 10203 5.02 48.3 

6 3 IIi 9017 4.16 71.6 8870 4.19 69.9 9049 4.21 68.8 



within 1-2%. The arithmetic mean of the empir- 
ical diagonal SO matrix elements (A s ° + A s _?)/2 
(see Fig. |5Jt) agrees with the ab initio SO splitting 
A so from Ref. [TJ] within few cm 1 near a mini- 
mum of these functions; (see Eq. ([8]). However, the 
present empirical A±(r) functions evaluated near 
the equilibrium distance r e of the triplet &-state 
are smaller on 20-25 cm -1 than their ab initio 
counterparts corresponding to the pure (c) Hund's 
coupling case [l|| (see TablellVf . Near the crossing 
point of singlet and triplet states r c sa 5.1 A(Fig. 
[7| the empirical off-diagonal SO coupling function 
£Xf>o coincide with ab initio result [14[ within 1 — 2 
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cm" 1 (Fig.lt)). 

The empirical matrix element -BCam (see Fig. 
[9^) reasonably agree with the ab initio estimate 
obtained in the framework of the second order per- 
turbation theory 



BL 



A-B^B-b 



Abl 



U B 



(7) 



[U A + U h )/2 

by using of the ab initio PEC for the B 1 Ii state 
Ub as well as SO £r° ?, and angular L\_ B coupling 



matrix elements [14j . Furthermore, the empirical 
non-equidistant SO splitting of the 6 3 II-statc (Fig. 
I and [5]) qualitatively agree with the estimates 



A* 



A so ± 5 A 3 



5A S 



E 



U b - U i 



(8) 



where the required PECs Uj for the j g 
a 3 E + ; c 3 E + ; B 1 !! states and relevant SO matrix el- 
ements A so , £1° of the 6 3 II state were taken from 
Ref. Hi. 



B. Details of the A E + -state interaction with 

the & 3 Iln=o,i sub-states 

The magnitude of the off-diagonal homogenous 
spin-orbit A 1 Y I + ~ 6 3 iIf2 = o interaction, see Fig. 
1 is comparable with vibrational spacing of the 
interacting states, see Table HVl This leads to sig- 
nificant regular non-adiabatic shifts of the most 
experimental level positions of the A ~ b complex 
with respect to their deperturbed (diabatic) sin- 
glet and triplet counterparts. Fig. ITUk gives us 
a typical example of the experimental level posi- 
tions of the complex given in a reduced scale as 
function of rotational quantum number J'. As 
can be seen, the experimental J'-dcpcndence of 
A state vibrational level identified for small J' as 
v* A = 1 (depicted by graph III), is shifted towards 
higher energies by approximately 75 cm -1 , or by 
about l.bu^ (compare with diabatic va = 1)- The 
graphs I - II belong to experimental = 20 and 
21, respectively. The experimental levels are sig- 
nificantly shifted down with respect their dibatic 
(deperturbed) positions. All three graphs I - III 
demonstrate a different slope from any of diabatic 
J'-dependencies. Moreover, the slope is chang- 
ing with J' demonstrating the singlet-triplet anti- 
crossing effect which takes place between experi- 
mental terms II and III around J' = 80 while be- 
tween I and II around J' = 130. 
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FIG. 7: (Color online) (a) The present empirical (em 
and available ab initio (abl-[lj^; ab2-\l^; a63-[lj|) 
PECs of the t4 1 £ + and &' ! iln=i states corresponding 
to the a Hund's coupling case. The insert zooms in the 
respective PECs around r e . (b) The present empirical 
diabatic PECs of the 6 3 IIn = o,i,2 sub-states and 

the corresponding adiabatic PECs (2; 3)fi = + in the 
c coupling case states. 



Along with the dominating A 1 T 1 + ~ £> 3 IIn=o in- 
teraction, weak local heterogenous perturbations 
caused by the (A — b)n=o ~ b 3 Hn=i indirect in- 
teraction are pronounced as sharp irregularities of 
the order of 1 cm -1 in the vicinity of intersection of 
experimental term values with diabatic dependen- 
cies of the 6 3 ITq = i sub-state. One such example is 
given by the inset in Fig. UOh while the more de- 
tailed analysis of the heterogenous perturbation is 
shown in Fig.[TlJ The peculiarity of the both cases 
is connected with experimental observation of two 
pairs of lines with the same J'-value, namely, for 
J' = 113/114 and J' = 157/158 (see Fig.Eb for the 
second example), giving levels with slightly diffcr- 
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FIG. 8: (Color online) Comparison of the present em- 
pirical first order spin-orbit functions with their ab ini- 
tio counterparts (ab) from Ref. [JJ]]. (a) the diago- 
nal empirical A±(r) and ab initio A 3 °(r) SO splitting 
functions; (b) the relevant off-diagonal f22,o( r ) SO cou- 
pling functions. r c is the internuclear distance where 
diabatic PECs intersect. 



ent energies. This is explained by a local pertur- 
bation, as a result of which, in each pair of lev- 
els with the same J', see Fig. [TTa and the inset 
in Fig. [TOb . one level is predominantly singlet A 
state, while another one is the triplet & 3 Iln=i sub- 
state. It should be noted that at intersection of 
graphs I and III with diabatic b\ graph in Fig. [TOa 
only small shifts were observed, without appear- 
ance of pairs with the same J'. This typical situ- 
ation should be attributed to too small fraction of 
singlet wavefunction in the partner level to observe 
it. 

Both types of interaction are described remark- 
ably well by the present deperturbation model, as 
is demonstrated by the residuals between the fit- 
ted and measured termvalues, see, for example, 
Fig.fTUb and Fig.lllb. Furthermore, the calculated 
fraction partitions Pi (see Fig.fTOb. Fig. [Tib. Table 
U and EPAPS [36]) unambiguously prove a domi- 
nant role of the direct SO A 1 Y> + ~ & 3 IIc> = o inter- 
action. Indeed, one can see that smooth variation 
of the singlet fraction in the total wave function 
corresponds to the smooth change of the respective 
termvalues. The sharp resonances at J' = 113/114 
(Fig. [TOt) and J' = 157/158 (Fig.QjV/) both corre- 
spond to the intersection of ./'-dependence of the 
respective experimental tervalues with diabatic de- 
pendence for Vbi — 17 caused by indirect interac- 
tion of A state with the 6 3 IIo = i sub-state. The 
very small irregularities of partition fractions no- 
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FIG. 9: Comparison of the present empirical and ab 
initio second order spin-orbit functions, (a) the em- 
pirical BC,Abi entering Eq. (0 while ab initio function 
is determined by Eq. ([7]); (b) the empirical SO split- 
ting functions A± entering Eq. Q while ab initio A± 
values are calculated by Eq. (JHJ . 



ticed on Fig. [TOb at J' = 81 and 124 indicate the 
predicted influence of 6 3 Il2 component which has 
only minor effect on energy and, hence, it could 
not be observed experimentally. 



C. Term values of the 41 K 133 Cs isotopomer 

To confirm mass-invariant properties of the de- 
perturbed parameters the rovibronic termvalues of 
the 41 K 133 Cs isotopomer have been predicted and 
compared with the experimentally observed ones. 
In the calculation only corresponding reduced mass 
fi of the isotopomer in the operator kinetic energy 
Eq. UJ and parameter B in the Eq. (J4| of the po- 
tential energy matrix defined by Eq. ([2|) and j3]) 
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FIG. 10: (a) Fragment of J'-dependence of the exper- 
imental term values as well as their deperturbed dia- 
batic singlet and triplet counterparts. Termvalues are 
given by the reduced scale: AE = E- 0.0222 x J'(J' + 
1). The graphs I, II, and III denote experimental data. 
Almost horizontal solid lines refer to diabatic singlet 
A S + state in the region va = 1 — 3; the dashed and 
dotted lines refer to diabatic triplet sub-states 6 3 IIn = o 
(v b0 = 18 - 20) and 6 3 ILn = i (v bl = 17 and 18), respec- 
tively. For small J' the effective vibrational quantum 
numbers can be attributed as v bQ (l) = 20, v b0 (II) = 
21 and v\(lll) = 1. The inset zooms in the splitting 
caused by the interaction with Vbi = 17 term. Ar- 
rows mark places of observed weak local perturbations 
with 61 sub-state. The triangles mark the levels which 
give the LIF progressions analyzed in Fig. \T3\ and 1141 
(b) Residuals of the fit Ej Xpt - Ef° for all respective 
J' values, (c) Fraction partition of the relevant non- 
adiabatic wavefunctions Pi for the experimental term 
values presented by graph II. The sharp resonance at 
J' = 113/114 corresponds to the intersection shown 
by inset in (a). The arrows in (c) mark the predicted 
influence of i> 3 Iln=2 component. 



was substituted. The calculated termvalues coin- 



FIG. 11: (Color online) (a) The experimental ter- 
mvalues around the local J' = 157/158 perturba- 
tion (see Fig. [5p) caused by the Vbi = 17 level of 
the diabatic & 3 IIn = i sub-state (marked by the solid 
line). Termvalues are presented by the reduced scale: 
AE = E - 10266 - 0.0237 J' (J' + 1). (b) Residuals of 
the fit E expt - E cc for respective J' values, (c) The 
respective partition functions Pi. 

cide with the experimental data (see Table [I]) with 
standard deviation of 0.0055 cm -1 . In absence of 
any adjustment of the parameters, the agreement 
even for the highly excited and strongly mixing 
levels can be assumed as excellent. 



D. Intensity distribution in the 

A 1 ^ ~ fe 3 n -> X 1 ^ LIF progressions 

As well-known a nodal structure of non- 
adiabatic vibrational wavefunctions is rather sen- 
sitive to strong intramolecular perturbations [l6l . 
l2ll . |22| . Therefore, a comparison of the experi- 
mental relative intensity distributions in a band 
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structure A ~ b — > X(ujf) LIF progressions with 
their theoretical counterparts would be a critical 
independent test on a completeness of the above 
energy-based deperturbation analysis. 

The theoretical transition probabilities lA~b^x 
were estimated as: 



■b->X 



(9) 



i/A-6-A- = £ CC (J') - Ex{v x ; Jx = J' ± 1) 



•b-X 



\dAx\vx) 



where dAx(f) is the ab initio A — X transition 
dipole moment from Ref's [H, EH- The power 
of 4 used in the transition wavenumber va~i>-*X 
assumes that the detector is proportional to the 
intensity of the incoming fluorescence light. The 
rovibronic energies Ex and eigenfunctions \vx) of 
the ground X 1 ^ state were obtained by solving 
the single channel radial equation with the em- 
pirical PEC from Ref. PH. Hereafter, both ex- 
perimental and calculated intensities obtained for 
P and i?-branches separately were then averaged 
since they basically were very close to each other. 

The remarkable agreement observed (see Fig. 
112a ) between the simulated and experimental in- 
tensity distribution in the very long vx €E [0, 80] 
LIF progression originating from the high rovi- 
bronic level of the complex convincingly proves a 
high reliability of the deperturbation model used. 
It is easily seen from Fig. [T2b that the non- 
adiabatic wavefunction <j> A (r) of the singlet A- 
state, which is responsible for the complicated os- 
cillation behavior of intensity observed on Fig.[T2"a. 
is only weekly perturbed by the triplet 6-state. The 
respective fraction of the singlet state Pa is domi- 
nating here being 0.81. 

To elucidate influence of the strong off diagonal 
homogenous SO A l H + ~ 6 3 IIn=o interactions on a 
nodal structure of the perturbed wavefunctions we 
have measured and simulated respective intensity 
distributions in A ~ b — > X{vx) LIF progressions 
originating from two rovibronic levels (marked by 
solid triangles on Fig. [TUb ) which are close to each 
other by J 1 and term values. These levels have, 
however, a quite different admixture of the b 3 Ho 
component. The agreement between calculated 
and experimental intensities is found to be, once 
more, remarkably good (see Fig. [T3"a.b). For both 
levels the I(vx) functions (Fig. [T3b.b) obey very 
well to so-called "Condon reflection approxima- 
tion" [r| since the vx dependence of I(vx) mimic 
the r-dependence of the initial wavefunction of the 
complex, ^(r)! 2 ; see Fig. [T3c. Indeed, in accord- 
ing to three lobes of non-adiabatic wavefunctions 



4>a (r) of singlet fraction depicted on Fig. [T5b the 
respective intensity distributions show three pro- 
nounced maximums which are clearly seen on Fig. 
[T3"a,b. At the same time, due to the strong ho- 
mogenous interaction an amplitude of the lobes of 
two wavefunctions are well distinguished leading 
to quiet different intensity distribution functions 
(compare Fig. [TBk with Fig. [T5b). It is interest- 
ing that the most dramatic changes in a shape of 
non-adiabatic wavefunctions are observed near the 
crossing point r c of the interacting states (Fig. [7b). 
For instance, the 4>a(i") of the J' = 44 level has 
very small additional lobe near r ~ 5.2 A. This 
lobe provides small but not negligible peak located 
at vx = 16 on the intensity distribution function; 
see Fig. Iblb . 

It also should be noticed that in according to 
11 oscillation theorem" |37[ the J' = 49 level men- 
tioned above would be assigned to the Va = 2 vi- 
brational state since its wavefunction <j>A(f) has 
two nodes at least. However, from the energy 
viewpoints (Fig. \Wk ) and Pa analysis it can be 
assigned as ^ = 1. The detailed analysis of this 
effect (conventional oscillation theorem breakdown 
caused by strong homogeneous perturbations) has 
been done separately in our accompanied paper 

& 

To investigate the influence of the local heteroge- 
nous (A — 6)n = o ~ b 3 Tln=i interaction on a nodal 
structure of wavefunctions we have measured and 
simulated intensity distributions in the LIF pro- 
gressions originating from two mutually perturbed 
levels, shown on Fig. [TTk . with the same J' = 158 
values. The agreement between calculated and ex- 
perimental values, once more, is remarkably good 
(Fig. [T4"a.b). It is surprising to see that the in- 
tensity distributions for both levels are found to 
be almost identical; the one level (marked as "s") 
having the lower energy demonstrates a dominant 
singlet A 1 Yi + character while the higher-lying level 
(marked as "t") has significant 6 3 no=i component. 
The calculation of the respective non-adiabatic 
wavefunctions <$> s Jf (r) (see Fig. [T4b ) confirms that 
their nodal structure is coinciding while the am- 
plitudes are different. Indeed, the wavefunctions 
relates to each other as 
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FIG. 12: (a) The experimental and calculated intensity 
distribution in the long LIF progression (see Fig. [5| 
originating Irom the high rovibronic level EA~b(J' = 
62) of the A ~ b complex, (b) Respective fraction of 
the non-adiabatic <j>A(r) wavefunction of the A-state 
(Pa = (M<Pa) = 0.81). 



E. Simulation of the optical cycle 

a 3 E+ -» A X E+ ~ & 3 n -»■ X J E+ 



Since the predictive abilities of the derived de- 
perturbation parameters of the A ~ 6 complex 
were unambiguously proved above we have used 
them to simulate a wavenumbers 
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FIG. 13: The experimental and calculated inten- 
sity distributions in LIF progressions originating from 
close-lying rovibronic levels of the A ~ b complex 
(marked by the solid triangles in Fig. 1 10 p : (a) for the 
vl = 21; J' = 44 having an admixture of the triplet 
6 3 IIn = o component (P b0 = (<f> b o\<f>bo) = 0.52); (b) for 
the v\ = 1; J' =49 level having a dominant singlet 
A*E + character (Pa = 0.74). (c) Respective non- 
adiabatic wavefunctions 4>a(t) of the singlet A-state. 



for the PUMP-DUMP optical cycle a 3 S+ -> 
A 1 Yi + ~ & 3 I1 — » which has been proposed in 

Rcf. H to transform an ultracold K(3 2 S)+Cs(6 2 S) 
colliding pairs in their absolute ground state. The 
excitation a — > A ~ & and emission A ~ 6 — > X 
transitions presented in Table fVl correspond to the 
most pronounced transition probabilities of the 
two-steps conversion process. Table [V] also con- 
tains the predicted radiative lifetimes TA~b (3&| 
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E(J'=158) = 10863.060 cm (P A =0.23) 
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FIG. 14: The experimental and calculated intensity 
distributions in the LIF progressions originating from 
two mutually perturbed levels of the A ~ b complex 
with the same J' = 158- value (see Fig. Illb ) : (a) hav- 
ing dominant sing let A X E + character ( P% = 0.58 
P b % = 0.14; P b \ = 0.28 ), marked in Fig. [TTJi by "s" 

(b) having dominant & 3 IIn = i component (P\ = 0.23 
P b f = 0.06; P b \ = 0.71 ), marked by "t" in Fig. [TTJi- 

(c) Respective non-adiabatic s /' (r) wavefunctions of 
the A-state. 



The wavenumbers u^^a 



cm 1 ) and transition moments M^ b _ x , 



DUMP (■ 
6i v A~b^X \ m 



M%*EZ (in 



a. u.) predicted for the most favorable stimulated Ra- 
man process a 3 E+(iV a = 0) -> 4 X E+ ~ b 3 Tl(J' = 
1) -> X^+ivx = 0; = 0). The ^5^6 values are 
given with respect to the common dissociation energy 
of the ground singlet X and triplet a states (see Fig.l). 
The ta~6 (in ns) and 7?a~6 (in %) are the radiative 
lifetimes and branching ratios of the spontaneous emis- 
sion, respectively. 
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and branching ratios of the spontaneous emission 

RA~b 

87T 2 g 2 

from the intermediate level of the complex with 
J' = 1 to the absolute ground level. Here, AUij = 
Uj — Uj is the difference potential. Both TA~b and 
RA~b values are also useful for estimating the effi- 
ciency of the simple PUMP-spontaneous emission 
process for producing population in the absolute 
ground level. 

The probabilities of PUMP transition were es- 
timated using the ab initio dipole transition mo- 



ment dba (r) from Ref. [TJ] and adiabatic wavefunc- 
tion \v a ) of the initial a 3 E + state [ll| calculated 
for a virtual vibrational level with the zero bound- 
ing energy E a (N = 0) = 0. The singlet-triplet 
X 1 T 1 + ~ a 3 E + perturbation of the initial level 
caused by a hyperfine Fermi contact interaction 
[TT| was neglected. The upper limit of the ter- 
mvalues for the A ~ b complex E^ C {J' = 1) was 
limited by 13200 cm" 1 in order to avoid an ex- 
trapolation error outside the experimental region. 
The wavenumbers and transition probabilities for 
arbitrary rovibronic levels of the A ~ b complex 
and ground singlet X and triplet a states can be 
generated by a request. 
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VI. CONCLUSIONS 

We have accomplished both experimental and 
deperturbation studies of the fully mixed A 1 Yi + 
and 6 3 II states of the KCs molecule based on a di- 
rect reduction of highly accurate Fourier-transform 
spectroscopy data on rovibronic termvalues of the 
A ~ b complex to potential energy curves and 
spin-orbit coupling functions for the mutually per- 
turbed states. 

The experimental data field starts from the 
lowest vibrational level v* A — of the singlet 
and nonuniformly covers rotational quantum num- 
bers J £ [7, 225] in the energy range E J £ 
[10040,13250] cm" 1 . Overall 42 adjusted fitting 
parameters have been required to reproduce more 
than 3400 rovibronic termvalues of the A ~ b com- 
plex with a standard deviation of 0.004 cm -1 which 
is absolutely consistent with the uncertainty of the 
experiment of 0.003-0.01 cm" 1 . Reliability of the 
derived structure parameters was unambiguously 
confirmed by a good agreement of the predicted 
termvalues of 41 K 133 Cs isotopomer and relative in- 
tensity distributions in the A ~ b — ► X(vx) LIF 
progressions with their experimental counterparts. 

Besides of the dominating homogeneous spin- 
orbit J 4 1 E + ~ 6 3 IIo = o interactions the local het- 
erogenous A 1 'E + r~j 6 3 ilf2=i perturbations have 
been discovered in the LIF spectra and their im- 
pact on a nodal structure of the non-adiabatic vi- 
brational wavefunctions of the complex have been 
analyzed by means of specially measured intensity 
distributions. It was found that the homogeneous 
perturbations can change dramatically a shape of 
wavefunctions (including even total numbers and 
position of their nodes) whereas the heterogenous 
perturbations affect only on amplitude of the wave- 
function in according to its fraction partition. 

The systematic deviations greater then 0.01 
cm" 1 of the experimental and reproduced termval- 
ues are still observed mostly for lying higher then 
12000 cm" 1 rovibronic levels of the complex es- 



pecially when they poses the significant admix- 
ture of the 6 3 IIn=i component. This drawback 
of the deperturbation model could be attributed 
to the neglected influence of first order spin-orbit 
and electronic-rotational interactions of the A ~ b 
complex with nearest c 3 E + state (Fig. [J). 

The strong A 1 S + ~ b 3 H spin-orbit interac- 
tion opens a window for direct observation of the 
formally spin-forbidden and basically very weak 
6 3 II — A" 1 £ + transitions originating from the rovi- 
bronic levels of the 6-state lying below va = 
of the perturbing A 1 T I + state. This work is in a 
progress. 

The Table fVl shows that the efficiency of the pro- 
posed conversion cycle a 3 T, + (N = 0) — > A 1 S + ~ 
6 3 n(J' = 1) -> X l H+{v x = 0; J x = 0) is lim- 
ited by its absorbtion part which requires a rather 
high power laser source generating in not very fa- 
vorable spectral diapason. This limitation caused 
by a very weak b — a triplet transition dipole mo- 
ment dba(r) in the considered r-range [H, [HI, has 
been discovered in the analogue cycle of NaRb [2l[ 
and NaCs [22j molecules. Nevertheless, relative 
high efficiency of the single step PUMP- sponta- 
neous emission transitions observed for the partic- 
ular levels of the A ~ b complex is encouraging for 
planning such experiments in the future. 
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